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propellant  ignition  by  the  combined  action  of  exti _ 

reaction  between  a  solid  fuel  and  gaseous  oxidizer  with-] 
ignition  as  the  limiting  case.  Several  commonly  used  deumnans  ox 
ignition  are  examined  and  compared  in  the  light  of  the  surface  tempera¬ 
ture  history.  Conditions  are  established  under  which  the  definition 
ignition  can  be  expected  to  influence  ignition  time;  major  differences 
appear  at  high  heating  rates  between  the  ignition  times  based  on  con¬ 
stant  ignition  temperature  and  on  other  criteria.  The  effect  of  oxi¬ 
dizer  concentration  on  ignition  time  is  investigated;  at  low  pressures, 
the  manner  of  variation  of  concentration  is  unimportant,  while  at  high 
pressures,  variations  of  mass  fraction  and  pressure  have  different 
effects.  Finally,  it  was  found  that  chemical  kinetics  is  the  principal 
factor  in  determining  ignition  time,  while  mass  diffusion  rate  imposes 
a. limit  on  the  attainment  of  ignition  without  significantly  affecting 
ignition  time. 
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NOMENCLATURE 

2 

Chemical  rate  coefficient  (gm/cn  sec) 

Specific  heat  (cal/gn°K) 

2 

Coefficient  of  mass  diffusion  (cm  /sec) 

Activation  energy  of  surface  reaction  (cal/mol) 
Constants  (Eq.  29,  30) 

Molecular  weight 

Heat  evolved  in  surface  reaction  (cal/gm  of  oxidizer 
reacted) 

Universal  gas  constant  (cal/mol°K) 

Temperature 

Surface  temperature  at  instant  of  contact  between  gas 
and  solid  phases  (Eq.  27,  28) 

Mass  fraction  of  oxidizer 

Empirical  constant  (Eq.  21) 

Empirical  constant  (Eq.  2*4 ) 

Specific  enthalpy  (cal/gm°K) 

Order  of  surface  reaction  with  respect  to  oxidizer 
2 

Pressure  (dyne/cm  ) 

2 

Flux  density  of  external  energy  (cal/cm  sec) 

Time  (sec) 

Distance  (cm) 
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a  Diffusion  parameter  (a=rT0/piD**QY0) ,  d5jnensionless 

6  Activation  energy  parameter  (8=E/HT0),  dimensionless 

y  External  heating  parameter  [y=qexp( 8)/QBY0n] , 
dimensionless 

6  Radiant  absorption  parameter  [6=YrX2viT0/r2q(Y+l)  J  * 
dimensionless 

0  Dimensionless  surface  temperature  (0=T  -T  /T  ) 

s  o  o 

X  Coefficient  of  thermal  conductivity  (cal/cra  sec°K) 

Vi  Radiant  absorption  coefficient  (cm*1) 
p  Density  (gm/cm^) 

x  Dimensionless  tine  {t=tt  (q+QBY0nexp(-E/RT0)  ]2t/r2T02) 

Tx  Gas  phase  thermal,  responsivity  IT1*(X|P  C  )  J 
cal/cn2  deg  K  sec^ 

k 

T0  Solid  phase  thermal  responsivity  [ro=(X2P2C2)  ] 


?rma. 

cal/cm2  deg  K  sec5 
T  Total  responsivity  (r=ri+r2) 
SUBSCRIPTS 

1  Gas  phase  (except  as  noted  in  text) 

2  Solid  phase  (except  as  noted  in  text) 
*  Ignition  conditions 

s  Surface 
c  Cutoff 
o  Initial 
«-  Far  Prom  surface 
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1.  INTRODUCTION 


Heterogeneous  reactions  have  been  proposed  as  the  controlling  mecha¬ 
nism  in  the  ignition  of  solid  propellants  (Ref.  l).  In  Ref.  2,  an  ex¬ 
tensive  review  is  given  of  the  status  of  ignition  theory  in  which  it  is 
noted  that,  in  the  case  of  heterogeneous  theory,  there  were  at  that  time 
(1966)  no  published  results  of  numerical  solutions  of  the  equations  re¬ 
presenting  heterogeneous  ignition,  either  with  or  without  external  heat¬ 
ing.  Williams  (Ref.  3)  later  published  an  analysis  accounting  for  ex¬ 
ternal  heating;  in  effect,  his  work  could  be  considered  a  study  oi  hy^er- 
golic  ignition.  The  purpose  of  this  paper  i u  to  extend  the  preceding 
analysis  (l)  by  including  the  effect  of  subsurface  absorption  of  exter¬ 
nal  radiant  flux,  (2)  by  including  the  effect  of  unequal  initial  tempera¬ 
tures  of  the  solid  and  gas  phases,  and  (3)  by  presenting  solutions  for 
heating  of  an  opaque  solid  propellant  by  a  combination  of  surface  chemi¬ 
cal  heating  and  a  continuous  or  interrupted  external  stimulus.  The  solu¬ 
tions  will  be  examined  in  the  light  of  several  arbitrary  definitions  of 
ignition  in  order  to  establish  the  conditions  under  which  the  choice  of 
ignition  criterion  may  or  may  not  be  important. 


2.  MATHEMATICAL  MODEL 


We  consider  the  transient,  one-dimens icnal  system  of  Fig.  1,  in  which 
the  solid  propellant  occupies  the  half-plane  x<0  and  the  gaseous  oxidizer 
occupies  the  half-plane  x>0.  Initially,  the  propellant  and  oxidizer  are 
at  the  uniform  initial  temperatures  T^  and  T^,  respectively.  At  time 
zero,  the  gasuous  oxidizer  of  uniform  mass  fraction  YQ  is  brought  into 
contact  with  the  propellant.  Simultaneously,  a  flux  of  constant  magni¬ 
tude  q  cal/cm2  sec  is  turned  on  and  illuminates  the  Interface  continu¬ 
ously  or  for  a  definite  period  of  time,  being  absorbed  by  the  solid  ac¬ 
cording  to  an  exponential  decay  law.  The  temperatures  at  a  great  dis¬ 
tance  into  the  solid  and  gas  are  maintained  at  and  T^  respectively, 
while  the  oxidizer  mass  fraction  at  a  great  distance  from  the  interface 
is  maintained  at  Y  .  In  addition  to  the  external  heating,  chemical  heat¬ 
ing  is  generated  a?  the  interface  according  to  a  one-step  irreversible 
reaction  governed  by  an  Arrhenius  dependence  on  surface  temperature. 
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FIC.  1.  Schematic  Representation  of 
Ignition  with  Surface  Reaction  and 
External  FTjlx. 

At  t>0,  the  processes  which  occur  include:  (l)  heat  exchange  between 
the  interface  and  the  solid  and  gas  phases,  (2)  diffusion  of  oxidizer  to 
the  surface,  (3)  surface  chemical  reaction  between  gaseous  oxidizer  and 
solid  fuel,  and  (4)  diffusion  of  products  from  the  surface.  The  follow- 


Ing  governing  equations  for  this  system  are  derived  from  general 
servation  laws  and  constitutive  relationships  (see,  e.g.  Ref.  4): 

con- 

Gas  Phase  (x>0) 

Energy:  p^(3h^/3t)  = 

afdj/C^)  (3h1/3x)]/3x 

(1) 

Mass:  p^SY/St) 

3[P;lD(3Y/3x)J/3x 

(2) 

Solid  PhAic  fx<0) 

Energy:  pg(3h2/3t)  = 

3 [ ( Xg/Cg ) ( 3hg/3x ) J /3x  +  yq  exp(px) 

(3) 

Boundary  (x*0) 

Temperature:  =  Tg  or 

hl/Cl  =  h2/C2  +  T2-*T1- 

(M 

Energy:  (X2/Cg)(3h2/3x) 

-  (X1/C1)(3h1/3X)=  Oo^OY/ax) 

(5) 

Mass:  p,D  (3Y/3x)  *  BY11  exp(-E/RT  ) 

JL  S 

(6) 

Initial  (t*»0) 

h^O;  h2»0;  Y*Yq 

(7) 

Asymptotic 
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F(X)=exp(62  X /it  ) er f c ( 6 /A/V )  -I  (12) 

For  convenience  in  considering  the  cases  of  interrupted  heating,  the 
roles  played  by  e  and  t  were  interchanged  from  Ref.  3  so  that  now  x  is 
the  independent  variable.  The  parameter  a  is  a  measure  of  the  relative 
importance  of  heat  and  mass  diffusions  with  a  value  of  zero  representing 
infinite  mass  diffusion;  B  is  the  dimensionless  activation  energy;  y  is 
the  ratio  of  external  heating  rate  to  initial  chemical  heating  rate;  and 
6  is  proportional  to  p,  the  absorption  coefficient. 

Although  a  computer  program  was  written  to  solve  Eq.  9—12,  cnly  the 
case  of  total  surface  absorption  (6=®)  has  been  considered  in  the  pa¬ 
rameter  study.  The  terms  involving  6  vanish  for  this  case.  In  the  cases 
of  interrupted  heating,  xc  is  the  dimensionless  cutoff  time.  For  continu¬ 
ous  heating,  tc  is  set  equal  to  t  so  that  only  the  first  integral  in  Eq. 

9  remains.  Complete  definitions  of  the  mathematical  symbols  ore  given  in 
the  nomenclature. 

In  Eq.  10,  the  expression  in  curly  brackets  represents  the  total 
heating  rate,  y  being  the  applied  external  rate  with  the  remaining  terns 
accounting  for  chemical  heating.  The  entire  function  H(0,X)  represents 
chemical  heating  after  shutoff  of  external  flux. 


3.  QUALITATIVE  BEHAVIOR 

The  general  nature  of  the  surface  temperature  history  depends  upon 
the  relative  importance  at  each  instant  in  time  of  the  constant  applied 
external  flux  and  the  variable  surface  chemical  heating.  Figure  2  illus¬ 
trates  all  the  characteristics  which  may  occur  in  the  temperature- time 
curve.  The  points  indicated  by  "1”  are  first  inflection  points.  When 
yc0,  this  point  Is  a  manifestation  of  the  increasing  reaction  rate  with 
temperature  before  diffusional  effects  limit  the  concentration.  When 
Y>0,  the  first  inflection  represents  the  occurrence  of  significant  chemi¬ 
cal  reaction  relative  to  external  heating.  The  point  marked  "y"  repre¬ 
sents  the  instant  when  chemical  and  external  heating  rates  are  equal. 

(it  is  absent  when  y=0  and  occurs  at  t=0  when  y“l)-  The  temperature  under 
these  conditions  is  readily  calculated  when  a=0  by  equating  the  chemical 
heating  rate  to  the  external  heating  rate  in  Eq.  10.  Then 

ey  =  lnY/(e-lnY)  (13) 

From  Eq.  13,  it  is  seen  that  there  are  combinations  of  heatirg  rate  and 
activation  energy  for  which  the  chemical  heating  cannot  become  dominant, 
viz.,  0<lnY-  The  underlying  basis  for  this  behavior  is  the  asymptotic 
form  of  the  Arrhenius  function  so  that  the  limitation  indicated  by  Eq.  13, 
while  relevant  only  for  quite  low  activation  energies  (ca.  4  kilo-calories 
ner  mole)  when  a=0,  becomes  of  increasing  importance  for  low  diffusion 
rates  (high  a). 
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FIG.  2.  Schematic  Heating  Curve  for  Surface 
Temperature  vith  and  without  External  Flux. 

The  points  marked  "2"  in  Fig.  2  are  second  inflection  points  where 
the  effect  of  decreasing  concentration  on  chemical  leaction  rate  begins 
to  be  important.  Second  inflection  points  occur,  albeit  at  unrealisti¬ 
cally  high  temperatures,  even  when  a»0  because  of  the  asymptotic  nature 
of  the  Arrhenius  function.  The  opposing  effects  of  increasing  tempera¬ 
ture  and  decreasing  concentration  ultimately  lead  to  a  maximum  in  the 
reaction  rate  (except  when  a=*0)  which  is  denoted  by  "m”  and  follows  the 
second  inflection  point.  The  remaining  characteristic  is  the  asymptotic 
value  of  surface  temperature  for  the  cases  y=0  ard  a^O.  This  behavior 
has  been  noted  and  explained  in  Ref.  2  and  3  and  represents  a  balance 
between  surface  reaction  rate  and  mass  diffusion  rate. 


3.1  HYPERGOLIC  CASES  (y=0) 

Figure  3  depicts  schematically  the  effect  of  a  and  6  on  temperature¬ 
time  curves  for  bypergolic  heating  (y«0).  The  tine  scale  is  different 
but  uniform  for  each  curve  in  order  to  illustrate  the  effect  of  a  and  B 
on  the  shape  of  the  0  vs  t  curve  for  a  wide  range  of  parameters.  High 
activation  energies  (large  0)  and/or  high  mass  diffusivities  are  typified 
by  long  induction  periods  followed  by  a  rapid  rise  in  temperature  (Curves 
A-C).  As  the  activation  energy  decreases,  the  temperature- time  curve 
gradually  loses  this  two-stage  characteristic  (curves  D  and  E);  the  in¬ 
flection  points  vanish  (curves  F-H);  and,  in  the  limiting  case  of  zero 
activation  energy,  the  dimensionless  surface  temperature  rises  smoothly 
to  a  value  of  1/a.  Lover  mass  diffusivities  promote  this  observed  tran¬ 
sition  at  higher  activation  energies. 


3.2  CONTINUOUS  EXTERNAL  HEATING  (y*0) 

The  effect  of  superimposing  a  constant  external  flux  upon  the  variable 
chemical  flux  can  be  argued  qualitatively  even  though  the  nonlinearity  of 
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FIG.  3.  Schematic  Surface  Temperature 
Histories  Shoving  Effect  of  Parameters 
for  Non-zero  Values  of  a. 


the  problem  precludes  any  true  mathematical  superposition.  As  the  inten¬ 
sity  of  external  heating  increases,  the  precipitous  nature  of  the  curves 
becomes  less  evident;  the  point  of  equal  heating  rates  is  no  longer 
reached;  the  inflection  points  disappear  and  only  the  point  of  maximum 
chemical  activity  remains,  shortly  before  depletion  of  the  oxidizer. 

This  transition  with  increasing  external  heating  is  represented  by  the 
series  of  curves  in  Fig.  3  except  for  the  upper  ends  of  the  curves, 
which  continue  to  increase  with  time  after  depletion  instead  of  approach¬ 
ing  the  asymptote  1 la  as  in  the  hypergolic  case. 

3-3  INTERRUPTED  EXTERNAL  HEATING  (y^O) 

Figure  U  illustrates  typical  surface  temperature  vs  time  curves  when 
the  external  flux  is  interrupted.  Curve  A  represents  shutoff  at  a  tine 
before  chemical  heating  has  become  appr  'iable,  either  because  of  (l) 
insufficient  heating  time,  (2)  low  diffusion  rates  and  consequent  early 
depletion,  or  (3)  low  activation  energy  and  high  external  heating  rate. 

In  curve  B,  the  heating  curve  has  already  deviated  from  that  for  an 
inert  solid  but  as  yet  insufficient  heat  is  being  generated  to  sustain 
the  reaction  after  shutoff.  Curve  C  typifies  shutoff  at  a  time  when 
chemical  heating  has  become  at  least  equal  to  external  heating.  Indeed 
a  usually  necessary  condition  for  the  occurrence  of  an  increasing  tempera¬ 
ture  soon  after  flux  shutoff  is  the  attainment  of  equal  chemical  and 
external  heating  rates.  The  requirement  is  unnecessary  when  the  activa¬ 
tion  energy  is  lov,  the  mass  diffusion  rate  is  high,  and  lny  is  near  its 
upper  limit  of  B,  since  under  these  conditions,  the  knee  of  the  Arrhenius 
curve  and  its  effect  on  rate  of  temperature  rise  occur  at  a  relatively 
lov  temperature.  The  requirement  is  insufficient  in  that  if  the  shutoff 
temperature  exceeds  1/a,  the  temperature  either  decays  slowly  toward  1/a 
(curve  D)  or,  when  the  overheating  is  severe,  drops  rapidly  and  recovers 
after  a  relatively  long  time  (curve  E). 


— 


I 


Ik* 


! 


f  ■ 


KWC  TP  fc86fc 


FIG.  5.  Solutions  for  Heterogeneous  Ignition  Model  with  Infinite 
Mass  Diffusion  Rate. 
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vhich  is  indicated  in  the  figure  for  lov  values  of  B-  As  previously 
remarked,  equal  chemical  and  external  heating  rates  occur  at  x=0  when 
Y*l.  Hence,  on  the  logarithmic  scale  used  in  Fig.  5»  all  curves  (solid) 
representing  this  aspect  of  the  temperature  history  are  asymptotic  to 
log  y=0.  The  dashed  lines  represent  the  tine  corresponding  to  the  second 
inflection  point  or  the  tine  at  vhich  the  rate  of  temperature  rise  vas 
so  rapid  that  the  computational  scheme  could  not  follow  it  without  speci¬ 
fication  of  an  inordinately  small  tine  increment.  The  two  sets  of  curves 
deviate  from  one  another  significantly  only  at  low  'external  heating  rates 
and  lov  activation  energies.  The  broken  curves  represent  equal  values  of 
0,  the  fractional  rise  in  surface  temperature  above  the  initial  interface 
temperature.  For  clarity,  only  a  few  values  of  0  are  shown.  At  high  ex¬ 
ternal  heating  rates  the  curves  approach  the  same  asymptotic  value  of  t 
for  all  activation  energies,  a  value  given  by  the  surface  temperature  of 
an  inert  system,  viz.,  t=(tt0)2/4,  as  obtained  from  Eq.  9-12  with  large  y- 
This  is  to  be  expected  because, with  high  heating  rates,  the  contribution 
due  to  chemical  reactions  is  no  longer  significant  and  the  system  appears 
inert.  As  the  external  heating  rate  decreases,  the  curves  representing 
different  constant  temperatures  for  a  given  activation  energy  all  coa¬ 
lesce  with  cne  another  and  with  the  curves  representing  precipitous  tem¬ 
perature  rise. 

Points  indicated  by  circles  represent  heating  times  for  what  are 
called  go-no-go  cases.  These  are  typified  by  curve  C  in  Fig.  4,  for 
which  the  surface  temperature  after  cutoff  of  external  flux  recovers 
(passes  through  a  minimum)  within  a  post  cutoff  time  interval  equal  to 
the  initial  heating  t'me.  Representative  calculations  reveal  that  for 
a=0  the  required  heating  time  for  go-no-go  behavior  lies  between  the 
points  designated  by  y  and  2  in  Fig.  2  (solid  curve).  Exceptions  occur 
as  y  approaches  its  upper  asymptotic  value  as  previously  noted. 


4.2  FINITE  DIFFUSION  RATES  (a>0) 

As  the  mass  diffusion  rate  decreases,  the  quantitative  results  of 
the  preceding  section  are  modified  in  two  ways.  First,  for  the  same 
values  of  8  and  y*  the  values  of  x  are  increased  by  an  amount  up  to  1C?, 
the  maximum  effect  occurring  at  low  external  heating  rates,  lev  activa¬ 
tion  energies,  and  lov  diffusion  rates.  Under  such  conditions  the  cut¬ 
off  time  has  less  effect  on  the  form  of  the  post  cutoff  temperature 
histories,  so  that  the  go-no-go  criterion  becomes  less  definite.  The 
effect  of  increasing  a  is  to  shift  upward  the  curves  representing  equal 
heating  rates  and  second  inflection  points  by  relatively  the  same  amount 
without  significantly  changing  the  spread  between  them.  Upper  asymptotic 
values  of  constant  temperature  curves  are  not  altered;  only  the  lover 
portions  are  affected  where  they  Join  the  curves  representing  second  in¬ 
flection  points. 
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Tbe  second  and  most  dramatic  effect  of  lov  mass  diffusion  rates  is 
the  manner  in  which  the  solid  curves  in  Fig.  5  are  truncated  so  that  only 
the  lover  sections  remain.  Qualitatively,  the  curves  for  a  given  value 
of  8  (activation  energy)  'end  vhen  t  corresponds  to  a  temperature  greater 
than  1/a.  Since  constant,  temperature  curves  axe  shown,  the  approximate 
upper  limit  of  x  is  readily  established  for  a  given  value  of  a.  The 
horizontal  portions  of  the  constant  temperature  curves  are  given  by 
Tlim*(if/2a)^.  Owing  to  the  relationship  between  go-no-go  behavior,  sec¬ 
ond  inflection  point,  and  point  of  equal  heating  rates,  this  second  limi¬ 
tation  greatly  reduces  the  range  of  parameters  over  which  go-no-go  be¬ 
havior  is  observed.*4 


U.3  HYPERGOLIC  CASES  (a=0) 


IJypergolic  heating  represents  not  only  the  limiting  situation  with 
no  external  heating,  but  also  is  the  condition  under  vhich  the  mass  dif¬ 
fusion  has  its  largest  influence  on  various  aspects  of  the  heating  curves. 

Owing  to  the  logarithmic  scale  employed  in  Fig.  5,  zero  external  heating 
cannot  be  represented.  These  cases  are  shown  graphically  in  Fig.  6  for 
diffusion  parameters  in  the  range  0  to  5  and  for  second  inflection 

points  ( T2 )  and  constant  temperature  (0  curves).  The  values  of  r  for  * 

8=0,  obtainable  analytically  for  a  given  6  (see  Ref.  2,  Eq.  A2*0  are 
used  as  an  aid  in  extrapolating  the  constant  d  curves  to  8=0.  In  the 
present  notation,  the  cited  equation  becomes  a6=l-exp(a2x/iT  )erfc(a2x/TT )  5. 

For  the  slow  mass  diffusion  rate  exemplified  by  a=10,  the  relationship  « 

of  surface  temperature  vs  tine  exhibits  no  second  inflection  point  except 
for  activation  energy  parameters  (8)  greater  than  50. 


5-  RELEVANCE  TO  PROPELLANT  IGNITION 

In  the  foregoing  sections  of  the  report,  attention  has  been  focused 
upon  formulating,  solving,  and  describing  a  model  in  vhich  heterogeneous 
(interface)  chemical  reactions,  augmented  by  an  external  interruptible 
energy  flux,  decemine  the  surface  temperature  history.  Granted  the 
adequacy  of  the  model,  one  is  still  faced  with  identification  of  a  given 
heating  curve  as  representing  ignition  or  nonignition  of  a  propellant. 
Mathematically,  the  choice  is  likely  to  be  highly  subjective,  especially 
for  situations  in  vhich  the  heating  curve  exhibits  no  outstanding  singu¬ 
larities.  Experimentally,  the  choice  nay  be  an  impossible  one  because  of 
the  difficulty  of  relating  experimental  observations  to  points  on  the 
heating  curve.  Convenient  experimental  procedures  include  detection  of 
light  and  pressure  rise  and  the  performance  of  go-no-go  tests. 


4  The  failure  to  attain  a  go-no-go  situation  is  an  artifact  of  the 
definition  of  go-no-go,  since,  if  the  surface  temperature  is  raised  above 
1/a  during  heating,  it  must  decrease  toward  1/a  (the  asymptotic  value) 
after  shutoff. 
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(a)  (b)  (c) 

FIG.  6.  Dependence  of  Dimensionless  Hypergolic  (y=0)  Ignition 
Time  upon  Activation  Energy  (8)  and  Diffusion  Hate  (a).  (a) 

a*0  and  0.1,  (b)  a=1.0,  (c)  a=5-0. 

We  shall  avoid  a  resolution  of  the  problem  of  defining  ignition  and 
instead  shall  examine  the  results  of  the  mathematical  solution  of  the 
model  in  the  light  of  several  commonly  used,  arbitrary  choices.  Speci¬ 
fically,  ve  shall  be  interested  in  constant  temperature,  equal  chemical 
and  external,  heating  rates,  second  inflection  point,  and  go-no-go.  We 
shall  further  wish  to  determine  the  conditions  under  which  the  defini¬ 
tion  of  ignition  is  of  minor  importance. 


5.1  CASES  WITH  EXTERNAL  HEATING 


For  convenience  in  the  ensuing  discussion  ve  designate  ignition 
times  based  on  ignition  temperature  as  t  ,  on  equal  chemical  and  exter¬ 
nal  heacing  rates  as  t  ,  on  second  inflection  point  as  t2*  and  on  inter¬ 
rupted  heating  (go-no-go)  as  Tj.  Reference  to  Fig.  5  reveals  the  exist¬ 
ence  of  two  regimes  of  external  heating  rate.  In  the  low  range,  charac¬ 
terized  by  y<100,  the  various  definitions  of  igniticn  lead  in  general  to 
different  ignition  times  according  to  the  inequalities 


<_  T2  (for  realistic  0) 


(15) 


It  is  noted  that  rY  <  T2  always,  and  that  Xj  is  closer  to  xY  at  low 
values  of  8.  The  relation  of  t q  in  this  scheme  depends  on  the  value  of 
0  chosen  for  the  ignition  criterion.  In  general,  x@  increases  with  0, 
but  the  dependence  becomes  very  weak  as  0  increases  and  y  decreases.  We 
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would  expect  this  trend  because  these  are  Just  the  conditions  favorable 
to  a  precipitous  rise  in  temperature.  The  effect  of  mass  diffusion  rate 

(a)  complicates  the  above  relations;  no  general  statements  can  be  made 
because  of  interactions  of  a,  8,  and  y. 

When  the  external  flux  corresponds  to  y^OO*  the  various  ignition 
times  become  less  distinct  so  that 

i.  t  ■  Tj  »  t2  (16) 

An  exception  to  Eq.  16  occurs  as  y  approaches  exp(B).  Then  tj  <  r  as 
indicated  in  Fig.  5  for  8=10.  When  y  exceeds  exp(8),  Ty  does  not  Y exist. 
The  regime  y>100  may  be  further  divided  accordingly  as:  (a)  Tg  =  x  ; 

(b)  Tg  <  T^;  (c)  tq  -  constant.  It  is  convenient  in  this  regard  to^refer 
to  two  special,  values  of  y.  These  are  Yj ,  below  which  x0  =  x  (constant 
temperature  lines  merge  with  equal  heating  rate  lines)  and  Y2  »  above 
which  tq  is  constant  ( const ant  temperature  lines  are  horizontal).  Values 
of  Yj  and  y2  are  shown  in  Fig.  7  and  8  as  functions  of  0  and  8  for  the 


FIG.  7*  Dimensionless  Relationship  Among  Surface  Temperature 
(0),  Activation  Energy  (8),  and  y\  (see  text)  for  Infinite 
Mass  Diffusion  Rate  (a=0). 
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FIG.  8.  Dimensionless  Relationship  Among  Surface  Temperature 
(8),  Activation  Energy  (0),  and  Y2  (see  text)  for  Infinite 
Mass  Diffusion  Rate  (a=0). 

case  a=0.  The  effect  of  a  is  twofold  as  was  notrd  in  connection  with 
Fig.  5*  First,  the  curves  of  y ^  (and  to  a  lesser  degree,  the  curves  of 
Yj )  are  shifted  slightly  to  the^left  as  a  increases.  Second,  the 
curves  are  truncated  so  that  only  the  sections  representing  0<l/a  remain. 
This  means  that  in  Fig.  7»  the  value  of  Yj  corresponding  to  given  values 
of  a(=l/9)  and  0  represents  the  highest  external  heating  rate  for  which 
go-no-go  ignition  is  observed. 

Between  y=100  and  Y=Yi »  all  the  definitions  of  ignition  considered 
here  lead  to  virtually  the  same  ignition  time.  This  range  is  narrow  for 
low  0  (activation  energies)  and  low-nass  diffusion  rates  (high  a),  van¬ 
ishing  completely  for  certain  combinations  of  a  and  0  as  inferred  from 
Fig.  7.  Under  such  conditions,  if  ignition  is  to  be  defined  at  all,  it 
must  be  based  on  a  critical  ignition  temperature,  ignoring  the  inevitable 
effects  of  quenching  by  thermal  and  diffusional  processes  if  the  stimula¬ 
ting  flux  is  terminated.  Above  y=Y2»  the  specification  of  a  critical 
ignition  temperature  is  tantamount  to  ignoring  all  chemical  and  diffu¬ 
sional  processes  and  treating  the  propellant  as  inert. 


5.2  CASES  WITH  NO  EXTERNAL  HEATING  (HYPERGOLIC) 

In  the  absence  of  external  flux,  the  times  and  Tj,  are  no  longer 
applicable  and  the  ignition  time  must  be  based  on  the  second  inflection 
point  or  a  critical  temperature.  Moreover,  there  is  now  a  limiting  value 
of  surface  temperature  rise,  given  by  8=l/a.  Hence,  for  low  diffusion 
rates  (high  a),  the  surface  temperature  is  limited  to  a  relatively  low 
value.  Even  the  adoption  of  a  critical  rate  of  rise  of  surface  tempera¬ 
ture  as  ignition  criterion  cannot  resolve  the  difficulty  in  this  instance 
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because  the  maximum  rate  occurs  at  zero  time.  The  relationships  between 
Tg  and  i2  are  shown  in  Fig,  6  for  various  combinations  of  a,  B*  and  0. 


0 


As  noted  at  ove,  values  of  0>l/a  are  not  attainable.  It  is  seen  that  x 
increases  with  0,  the  effect  being  greatest  for  small  values  of  B  and" 
large  values  of  a.  The  value  of  t2  may  be  greater  or  less  thai.  in 
contrast  to  cases  with  external  heating  because  the  second  inflection 
point  frequently  occurs  at  low  surface  temperatures  in  the  hypergolic 
case.  The  values  of  and  x2  approach  one  another  as  B  increases;  hence, 

hypergolic  ignition  may  be  equally  well  defined  by  some  critical  tempera¬ 
ture  (any  value  <l/a)  or  second  inflection  point  at  higher  activation 
energies  exceeding  the  value  Bi,  a  function  of  a  as  shown  in  Fig.  6. 

6.  INTERPRETATION  AND  IMPLICATIONS  OF  RESULTS 

It  has  been  shown  that  the  definition  of  ignition  influences  ignition 
time  depending  upon  the  range  of  parameters  involved.  Comparison  has 
been  made  in  terms  of  certain  dimensionless  parameters  and  variables 
which  were  convenient  in  the  numerical  solution  of  the  problem.  The 
results,  displayed  in  Fig.  5  and  6  may  be  indicated  functionally  by  the 
general  relation: 

t.  =  f  (a,  6,  y»  n,  Def)  (17) 

or  in  dimensional  form 

u[q+QBY  nexp(-E/RT  )]2t./r2T2  = 
o  o  o 

f[rTo/PjD,sQYo,  K/RTo,q  exp(E/RTo)/QBYon,  n,  Def]  (17a) 

where  Def  denotes  dependence  on  definition  of  ignition.  It  is  immediately 
obvfous  that  such  interesting  physical  quantities  as  q,  Y  ,  Q,  3,  E,  and 
D,  as  well  as  pressure  (which  enters  through  T,  p,  B,  D,  and,  under  cer¬ 
tain  conditions,  TQ)  are  involved  in  a  complicated  manner.  Moreover,  it 
is  not  possible,  by  the  choice  of  any  other  single  set  of  dimensionless 
groups,  to  decouple  all  the  physical  quantities  whose  effect  is  desired. 

In  this  section  the  effects  of  t^ss  fraction,  thermal  properties,  and 
external  heating  rate  will  be  aained  by  the  use  of  Fig.  5  and  6  and  the 
requisite  additional  graphical  relations.  The  discussion  will  be  in 
terms  of  the  magnitude  of  external  heating  and  the  ignition  criterion. 

The  more  complicated  effect  of  pressure  will  be  dealt  with  separately. 
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6.1  EFFECTS  OTHER  THAN  OF  PRESSURE 

6.1.1  Hyperbolic  Ignition  ( y=0) 

Under  conditions  of  no  external  heating,  only  two  meaningful  ignition 
criteria  remain  of  the  ones  previously  discussed,  viz.,  the  second  inflec¬ 
tion  point,  t2,  and  constant  temperature,  xQ.  Equation  17a  simplifies  to 

ir[QBYonexp(-E/RTo)/rTo]2t.  =  f  E/RTq,  Def]  (18) 

This  relationship  is  shown  graphically  in  Fig.  6  for  n*l.  The  effect 
of  n  is  second  order  for  a<5  and  8>20.  The  function  Mff*  in  Ec  •  18  is 
given  empirically  within  5Jf  by 

f  -  f  (E/RT  )  *  (19) 

JL  O 

where  f^  and  f^  are  two  functions  of  a  shown  below. 


0,  .1  2.3  2.23 

1  5-0  -2.UU 

5  6.9  -2.UU 

At  constant  E/RT  ,  tp  changes  by  less  than  a  factor  of  2  over  the  range 
of  values  of  a  investigated.  Ignoring  this  variation  compared  to  order 
of  magnitude  effects  allows  us  to  conclude  that 

t  “ (TT  / QBY  n ) 2 ( E/RT  ) 2 • 3  (20) 

d  o  o  o 

The  effect  of  a  would  be  to  increase  the  exponent  of  the  first  factor  of 
Eq.  20  very  slightly. 

If  ignition  is  defined  by  the  constant  temperature  criterion,  t  is 
given  by  Eq.  20  where  0>B,.  When  10<6<82»  for  meaningful  temperature 
rises  represented  by  0>O.2,  tQ  differs  by  no  more  than  30 i  from  the  value 
given  by  Eq.  20.  Under  no  conditions  can  ignition  occur  under  this  cri¬ 
terion  if  0.  >l/a. 

ign 

6.1.2  Ignition  with  External  Heating  (y>0) 

The  discussion  of  ignition  with  external  heating  will  be  limited  to 
the  practical  range  y>100,  thereby  serving  the  following  useful  purposes: 

(l)  Only  two  of  the  ignition' criteria  need  be  considered,  viz., 
constant  temperature  and  go- no-go. 
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Eq. 


(2)  The  chemical  heating  term  relative  to  q  on  the  left  side  of 
17a  is  negligible. 


(3)  The  effect  of  a  nay  be  ignored  except  as  it  establishes  limit¬ 
ing  conditions. 


(k)  The  direct  effect  of 
Y  )  may  be  ignored. 


in  solving  Eq.  9-12  (not  its  effect  in 


Figures  9 ,  10,  and  11  are  derived  from  Fig.  5  and  are  included  to 
show  the  individual  effects  of  several  physical  parameters.  Figure  9 
emphasizes  primarily  the  effect  of  external  heating  rate  on  ignition  time. 
In  order  to  clarify  the  presentation,  only  a  limited  number  of  constant 
temperature  curves  are  displayed.  The  trend  in  general  follows  the  one 
set  in  solid  phase  thermal  theory  in  which  higher  surface  temperatures 
are  required  for  ignition  at  higher  external  heating  rates.  Between 
Y*100  and  ysYj  there  is  little  to  distinguish  ignition  times  defined  by 
go-no-go  and  constant  temperature;  the  relationship  followed  is 


.-a 

q 


(21) 


where  ’’a”  is  a  weak  function  of  q  and  8,  and  is  given  over  a  wide  range 
of  parameters  by 

a  =  2-(8.  k/B)  (22) 

As  y  decreases  toward  100  or  increases  toward  exp  (8)»  "a"  is  signifi¬ 
cantly  lover  than  predicted  by  Eq.  22.  Above  y“Yj  *  go-no-go  ignition 
times  continue  to  follow  Eq.  21  subject  to  the  stated  limitations,  while 
ignition  times  based  on  constant  temperature  above  y=Y2  are  given  by 


t#  =  irr2(T,-T  )2Aq2 


(23) 


It  is  emphasized  that  Y]  Y2  depend  upon  8  and  0,  the  critical  tempera¬ 

ture  chosen  for  the  constant  temperature  criterion.  Equation  23  is  the 
classical  expression  for  ignition  tine  of  an  inert  solid  based  on  constant 
ignition  temperature  T#  (Ref.  2). 


Figure  10  depicts  the  effect  of  Q,  the  heat  of  reaction,  B,  the  rate 
coefficient,  and  Y  ,  the  initial  mass  fraction  of  oxidizer  upon  ignition 
time.  At  heating  rates  above  constant  temperature  criterion  pre 
diets  no  effect  of  the  group,  QBY0n,  while  the  relation  for  the  go-no-go 
criterion  becomes 


where 


(0BYonrb 


b  =  8.k/8 


(2k) 


(25) 
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FIC-  10.  Effect  of  QBYqi1  on  Ignition  Tine  (Constant  Tempera¬ 
ture  and  Go-IJo-Go  Criteria). 

inert  diluent  in  the  gas  phase.  It  is  also  seen  that  initial  surface 
temperature,  TQ,  has  a  powerful  effect  on  ignition  time,  especially  at 
high  values  of  QBY0n/q.  Figure  11  shows  the  importance  of  activation 
energy  in  determining  ignition  time. 

6.2  EFFECT  OF  PRESSURE 

Ambient  gas  pressure  influences  the  ignition  phenomenon  by  three 
basic  mechanisms; 

(l)  Effect  on  concentration  of  active  species  through  the  rate 
coefficient  B,  given  by 


B  =  vpjn  -  v(pM/RTs)n  (26) 

The  effect  of  Tg  on  B  has  been  neglected  relative  to  the  effect  of  Tg  in 
the  Arrhenius  factor  of  reaction  rate. 

(2)  Effect  on  thermal  losses  by  conduction  to  the  gas  as  reflected 
in  f  and  TQ. 
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FIG,  11.  Effect  of  Activation  Energy  on  Ignition  Time  (Constant 
Temperature  and  Go-No-Go  Criteria), 

(3)  Effect  on  diffusional  processes  in  the  parameter  a=rTo/p1DisQY0. 

Mechanism  (1)  is  of  importance  over  the  entire  pressure  range.  At 
low  pressures,  it  leads  to  the  conclusion  that  ignition  time  depends  upon 
oxidizer  concentration,  regardless  of  whether  it  is  varied  by  changing 
the  mass  fraction  or  t  :e  total  pressure. 


The  effect  of  mechanism  (2)  is  obtainable  from  the  expressions  for  T 
and  Tq: 

r  *  (XjPjCj)5*  +  (XjPjCj,)1*  (27) 

TT0  =  (XjPjCj  +  (X202c2)\m  (28) 

When  the  solid  and  gas  phases  are  initially  at  the  same  temperature, 

Tla>  -  T^;  the  initial  surface  temperature  is  then  independent  of  thermal 
properties.  If  Tleo  and  T2qo  are  unequal,  Tand  Tq  both  depend  on  thermal 
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properties;  hence  they  depend  on  pressure.  Incorporating  the  usual  rela¬ 
tions  {Ref.  U  )  between  pressure  and  gas  transport  properties  leads  to 


r  -  Kjp*5  +  k2 

(29) 

(k3p,*+k4)/(kip5s+k2) 

(30) 

The  quantitative  effect  of  pressure  on  ignition  time  through  the  heat  loss 
mechanism  is  obtainable  by  evaluation  of  Kj ,  K^,  K3,  and  K  for  a  particu¬ 
lar  situation,  followed  by  use  of  the  appropriate  graphical  or  empirical 
relationship  between  T,  TQ,  and  ignition  time. 

The  effect  of  pressure  on  diffusional  processes  appears  only  in  the 
parameter  a;  nence,  it  would  appear  only  in  a  limiting  role  under  the 
go-no-go  definition  of  ignition.  Qualitative  evaluation  of  this  effect 
of  pressure  depends  upon  establishing  the  pressure  dependence  of  each  of 
the  terms  in  the  definition  of  a,  as  well  as  the  pressure  dependence  of 
the  limiting  value  of  a  as  nay  be  inferred  from  Fig.  7.  The  complicated 
nature  of  the  relationships  involved  precludes  any  comprehensive  study 
or  general  conclusions. 

6.2.1  Hypergolic  Ignition 

The  first  order  effects  of  pressure  on  hypergolic  ignition  t5me  are 
obtainable  from  Eq.  20.  Two  cases  are  important,  viz.,  equal  or  unequal 
initial  temperatures  of  the  gas  and  solid. 

6.2. 1.1  Equal  Initial  Temperatures.  Here  the  value  of  TQ,  the 
initial  interface  temperature,  ir  pressure  independent,  so  that 

t.  -  (K1p?*+K2)2p-2n  Y0_2n  (31) 

At  low  pressures  it  is  clear  that  p  and  YQ  influence  ignition  time  in 
nearly  the  same  manner.  Thus  we  have  the  often  repeated  statement  that 
concentration  changes  affect  the  ignition  time  in  the  same  manner, 
whether  brought  about  by  changes  in  pressure  or  mass  fraction  of  the 
oxidizer.  At  higher  pressures,  the  effects  are  not  the  same:  the  sensi¬ 
tivity  of  t#  to  pressure  changes  is  less  owing  to  enhancement  of  heat 
loss  to  the  gas  phase  at  high  pressures.  For  an  effective  zero  order 
reaction,  there  would  be  no  effect  of  changing  the  mass  fraction,  while 
ignition  time  would  increase  with  increasing  pressure. 

6. 2. 1.2  Unequal  Initial  Temperatures .  Under  this  condition  we  may 
write  Eq.  20  as 


t.  -  (K1p,S+K2)2T0-*V2lV2n 


(32) 
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Again  the  effects  of  pressure  and  mass  fraction  are  similar  at  low  pres¬ 
sures.  At  elevated  pressures,  TQ  is  dependent  upon  pressure.  Two  situa¬ 
tions  may  be  considered  depending  on  whether  the  solid  or  gas  is  at  a 
higher  initial  temperature. 

(1)  Tlco  >  T2a>.  With  a  higher  temperature  in  the  gas,  TQ  increases 
with  increasing  pressure,  leading  to  a  greater  dependence  of  t#  on  pres¬ 
sure  than  was  found  in  Section  6. 2. 1.1  above. 


(2)  T2a>  >  Tc  decreases  with  increasing  pressure,  so  that  the 

effect  of  pressure  on  ignition  time  is  less  than  with  equal  initial 
temperatures . 

6.2.2  Ignition  with  External  Heating 

Under  conditions  involving  external  stimulus,  consideration  must  be 
given  to  the  definition  of  ignition  as  well  as  to  the  relative  tempera¬ 
tures  of  the  solid  and  gas. 

6.2. 2.1  Equal  Initial  Solid  and  Gas  Temperatures.  The  criterion  of 
an  ignition  temperature  (Lq.  23)  leads  to  the  result  that,  for 
changes  of  pressure  in  the  low  pressure  regime  do  not  affect  ignition 
time.  At  high  pressures,  owing  to  greater  heat  loss,  the  ignition  tine 
increases  with  pressure.  For  the  ignition  temperature  criterion  and 
y>YA,  or  for  go-no-go  and  y>100»  Fig.  10  may  be  applied  in  the  form  of 
the  relation 

t#  -  (KlP^K2)* 1 2 * * * * *(p/Yor8‘4nRTo/E  (33) 


At  low  pressures,  we  again  find  that  the  influence  of  pressure  and  mass 
fraction  on  ignition  time  are  similar,  although  the  sensitivity  is  less 
than  under  hypergolic  conditions.  At  high  pressures,  the  effect  of  in¬ 
creasing  pressure  could  conceivably  result  in  increasing  the  ignition 
time,  especially  for  reactions  with  high  activation  energies. 


6. 2.2. 2  Unequal  Initial  Solid  and  Gas  Temperatures.  Equations  23 
and  33  may  also  be  applied  for  the  temperature  and  go-no-go  criteria, 
respectively,  taking  into  account  the  dependence  of  TQ  on  pressure.  Two 
situations  arise,  depending  on  the  relative  magnitude  of  T1<jo  and  T2a>. 

We  consider  the  high  pressure  regime. 


(1)  T1<jb  >  T2o#.  The  initial  interface  temperature,  T  ,  increases 

with  increasing  pressure,  leading  to  decreasing  ignition  times  for  both 
go-no-go  and  temperature  criteria.  The  effect  of  pressure  on  the  go-no- 
go  ignition  time  is  extremely  great,  however,  owing  to  the  exponential 
nature  cf  TQ  in  Eq.  33*  -17 

(2)  T^  >  Tloo.  The  initial  temperature  decreases  with  increasing 
pressure.  Ir.e  ignition  time  based  on  constant  temperature  increases  with 
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increasing  pressure  while  the  effect  on  go-no-go  tine  is  indefinite  de¬ 
pending  on  the  relative  effects  of  the  two  factors  in  Eq.  33* 


7.  SUMMARY 


Results  have  been  presented  of  computations  on  a  mathematical  model 
of  solid  propellant  ignition  in  which  constant  external  heating  ard 
surface  chemical  reactions  provide  the  ignition  stinulus.  One  of  the 
reactants  (oxidizer)  is  assumed  to  be  initially  present  in  the  gas  phase 
and  to  react  with  the  solid  fuel.  Mass  and  heat  diffusion  and  conduction 
are  included  in  the  analysis  while  surface  regression  is  ignored  during 
the  ignition  transient.  The  surface  temperature  history  has  been  computed 
for  the  cases  of  continuous  and  interrupted  application  of  external  energy 
and  for  the  hypergolic  case  with  no  external  flux.  A  wide  range  of  values 
of  kinetic  constants,  diffusional  parameters,  and  heating  rates  was 
studied.  Several  commonly  used  definitions  of  ignition  have  been  examined 
in  terms  of  the  surface  temperature  history.  These  include:  (l)  attain¬ 
ment  of  a  specified  surface  temperature,  (2)  attainment  of  a  second  inflec¬ 
tion  point  in  the  curve  of  temperature  vs  time,  (3)  attainment  of  a  condi¬ 
tion  in  which  chemical  heating  equals  external  heating,  and  (h)  in  cases 
of  interrupted  heating,  the  occurrence  of  a  runaway  surface  temperature 
after  shutoff,  within  a  time  interval  equal  to  the  initial  heating  period. 
All  or  none  of  these  four  characteristics  of  the  heating  curve  ere  ob¬ 
served  depending  upon  the  set  of  physical  parameters  involved. 

Under  hypergolic  conditions,  only  the  first  two  definitions  are  rele¬ 
vant,  becoming  nearly  equivalent  for  moderate  tc  high  activation  energies 
and  mass  diffusion  rates.  For  low  mass  diffusion  rates,  the  maximum  tem¬ 
perature  rise  is  small,  so  that  even  the  appearance  of  an  inflection 
point  in  the  temperature  vs  tine  curve  is  no  assurance  of  physical  igni¬ 
tion.  For  low  to  moderately  high  external  heating  rates,  all  four  defi¬ 
nitions  lead  to  nearly  equal  ignition  times;  under  extremely  high  external 
rates,  the  first  definition  results  in  the  prediction  of  shorter  ignition 
times.  In  addition,  there  exists  an  upper  external  heating  rate  above 
which  the  fourth  definition  does  not  lead  to  ignition  at  all.  This  limit¬ 
ing  rate  is  an  increasing  function  of  mass  diffusion  rate  and  activation 
energy;  however,  when  ignition  occurs,  ignition  time  is  relatively  insen¬ 
sitive  to  mass  diffusion  rate. 

The  effect  of  physical  parameters  and  groups  of  parameters  has  been 
examined.  In  particular,  owing  to  conclusions  and  claims  of  previous  in¬ 
vestigators,  attention  has  baen  directed  toward  studying  the  effect  of 
oxidizer  concentration  on  ignition  time.  It  was  found  that  in  the  low 
pressure  range,  for  both  hypergolic  and  externally  heated  cases,  ignition 
time  is  nearly  the  same,  whether  concentration  is  varied  by  changing  the 
mass  fraction  of  oxidizer  or  the  total  pressure.  However,  at  higher  pres¬ 
sures  (exact  value  a  function  of  the  properties  of  materials  involved). 


22 


REFERENCES 


Anderson,  R. ,  R.  S.  Brown,  and  L.  J.  Shannon,  "ignition  Theory  of 
Solid  Propellants,"  AIAA  Preprint  6^-156  (196U). 

Price,  E.  W. ,  H.  H.  Bradley,  Jr.,  G.  L.  Dehority,  and  M.  M.  Ibiricu. 
"Theory  of  Ignition  of  Solid  Propellants,"  AKER  INST  AERONAUT 
ASTRONAUT  J,  Vol.  U,  No.  7  (July  1966),  pp.  1153-81. 

WiUians,  F.  A.  "Theory  of  Propellant  Ignition  by  Heterogeneous 
Reaction,"  AMER  INST  AERONAUT  ASTRONAUT  J,  Vol.  U,  No.  8  (August 
1966),  pp.  135U-57. 

Bird,  R.  B. ,  W.  E.  Stewart,  and  E.  N.  Lightfoot.  Transport  Phenomena. 
John  Wiley  and  Sons,  New  York,  i960. 


NWC  TP  H86U 


\ 

) 

1 

j 

< 


INITIAL  DISTRIBUTION 


11  Naval  Air  Systems  Command 
AIR-03  (1) 

AIR-03B  (1) 

AIR-310  (1) 

AIR-330  (1) 

AIR-330B  (1) 

AIR-503  (1) 

AIR-536  (1) 

AIR-5366  (1) 

AIR-5367  (1) 

AIR-601*  (2) 

1*  Chief  of  Naval  Material 
NSP-00  (1) 

NSP-001  (1) 

NSP-20  (1) 

NSP-2733  (1) 

5  Naval  Ordnance  Systems  Command 
0RD-033  (1) 

0RD-03311  (1) 

0RD-03312  (1) 

ORD-9132  (2) 


3  Office  of  Naval  Research 
Code  102  (1) 

Code  1*29  (1) 

Code  473  (1) 

2  Naval  Ordnance  Laboratory,  White  Oak 
Carl  Boyars  (l) 

Technical  Library  (l) 

2  Naval  Ordnance  Station,  Indian  Head 

L.  A.  Dickinson  (l) 

Technical  Library  (l) 

3  Army  Missile  Command,  Redstone  Arsenal  (Redstone  Scientific  Information 
Center,  Chief  Document  Section) 

2  Army  Ballistics  Resecrch  Laboratories,  Aberdeen  Proving  Ground 
Interior  Ballistics  Laboratory 
R.  C.  Strittmater  (l) 

Technical  Library  (l) 

2  Picatinny  Arsenal 
J.  Picard  (l)- 
Technical  Library  (l) 

1  Sunflower  Army  Ammunition  Plant  (SMUSU-R) 


25 


NWC  TP  U861* 


1  Air  Force  Office  of  Scientific  Research  (Propulsion  Division, 

Dr.  Bernard  T.  Wolf son) 

2  Air  Force  Rocket  Propulsion  Laboratory,  Edwards  (RPMCP)  ^ 

Capt.  Stamets  (l) 

Richard  Spann  (1) 

2  6593  Test  Group,  Edwards  Air  Force  Base  (Development,  Don  Hart) 

1  Wright-Patterson  Air  Force  Base  (APRT,  Lt.  R.  D.  Chronister) 

1  Advanced  Research  Projects  Agency  (Technical  Information  Office) 

20  Defense  Documentation  Center 

1  Bureau  of  Mines,  Pittsburgh  (Explosive  Research  Center,  I.  Liebman) 

19  National  Aeronautics  &  Space  Administration 
MGS,  E.  Hall  (l) 

MTA,  M.  C.  Waugh  (l) 

NaPO,  John  J.  Stashak  (l) 

RC,  J.  L.  Sloop  (l) 

RF,  R.  D.  Ginter  (l)  i. •• 

RP,  A.  0.  Tischler  (l) 

RPM,  William  Cohen  (1) 

RPS,  R.  Ziem  (5) 

RTA,  R.  V.  Hensley  (1) 

RTP,  J.  J.  Phillips  (1) 

RV,  M.  B.  Ames  (l) 

RV-1,  R.  W.  May,  Jr.  (l)  - 

SV,  V.  L.  Johnson  (l) 

ATSS-AC,  Technical  Library  (2) 

1  George  C.  Marshall  Space  Flight  Center  (Technical  Library) 

1  Goddard  Space  Flight  Center  (Technical  Library)  * 

2  Langley  Research  Center 

Robert  L.  Swain  (l) 

Technical  Library  (l) 

U  Lewis  Research  Center 

Dr.  Louis  A.  Povinelli  (l) 

Dr.  Richard  J.  Priem* (l) 

James  J.  Kramer  (l) 

Technical  Library  (l) 

2  Manned  Spacecraft  Center 

J.  G.  Thibodaux  (l) 

Technical  Library  (l) 

1  Aerojet-General  Corporation,  Azusa,  Calif.  (Technical  Library),  via  NPRO 

3  Aerojet-General  Corporation,  Sacramento,  via  AFPRO 

J.  Wiegand  (l) 

Department  U730  (l) 

Head,  Technical  Information  Office  (l) 

1  Aerospace  Corporation,  Los  Angeles  (Technical  Library) 

2  Allegany  Ballistics  Laboratory,  Cumberland,  Md.  ^ 

T.  A.  Angelus  (l) 

Technical  Library  (l) 

2  Applied  Physics  Laboratory,  JHU,  Silver  Spring 

Dr.  W.  it,  Avery  (l)  > 

Dr.  Robert  H.  Cantrell  (l) 


26 


£ 


;  i  ,  %  *>.■ 


r?  n 


-laaA^TFTFn 


Security  Classification 


DOCUMENT  CONTROL  DATA  -  RAD 


1.  ORIGINATING  ACTIVITY  (Cmtpormim  »<*•*> 

Haval  Weapons  Center 

China  Lake,  California  93555 

2m  RCRORT  »ICU«iTt  C  k  A»«<rtCA  TlOM 

U5ELASSTFTED 

tO  RAOtiP 

1  REPORT  TITLE 

Heterogeneous  Ignition  Theory 

4  DESCRIPTIVE  NOTES  (Tjrpm  • >1  rmpmrt  «rW  mcNtlf,  *#*•«; 

Pesearch  Report 

t  AUTHORS;  fU«l  namm.  #»r,f  nmmm  fntttmi) 

Bradley,  H.  H. ,  Jr. 

••  REPORT  OATE 

February  1970 

26  u 

|»  CONTRACT  or  «AAmT  mo. 

a  PROJECT  MO. 

-TA  0RD-331-001/20O-1  URO-100-202 

4t 

9 A.  ORI9MATORT  REPORT  MMMOERfSJ 

NVC  TP  1*86Ib 

9  A  OTMKR  R^AQRT  MOOj  (4f  .Aw  A»Wn  MMfW* A»I#mO 

IB  AVAIL  AR1LITY/LIMITATION  NOTICES 

This  document  is  subject  to  special  export  controls  and  each  transmittal  to 
foreign  governments  or  foreign  nationals  may  be  made  only  with  prior  approval 

II.  SUPPLEMEMTARV  motes 

12  SPORSORING  MILITARY  ACTIVITY 

Naval  Ordnance  Systems  Command 

Department  of  the  Navy 

Washington.  D.  C. 

tf.  ABSTRACT 

Mathematical  solutions  are  obtained  for  a  model  of  solid  propellant  Ignition 
by  the  combined  action  of  external  flux  and  surface  reaction  between  a  solid 
fuel  and  gaseous  oxidizer  with  hypergolic  ignition  as  the  limiting  case.  Several 
cocaonly  used  definitions  of  ignition  are  examined  and  compared  in  the  light  of 
the  surface  temperature  history.  Conditions  are  established  under  which  the 
definition  of  ignition  can  be  expected  to  influence  ignition  time;  major  differ¬ 
ences  appear  at  high  heating  rates  between  the  ignition  times  based  on  constant 
ignition  temperature  and  on  other  criteria.  The  effect  of  oxidizer  concentration 
on  *t longtime  is  investigated;  at  low  pressures,  the  manner  of  variation  of 
concentration  is  unimportant,  while  at  high  pressures,  variations  of  mass  fraction 
and  pressure  have  different  effects.  Finally,  it  was  found  that  chemical,  kinetics 
is  the  principal  factor  in  determining  ignition  time,  while  mass  diffusion  rate 
imposes  a  l*-ait  on  the  attainment  of  ignition  without  significantly  affecting 
ignition  time.^^  ^ 


DD  /J!-m  1473 


0101 *107-4200 


UNCLASSIFIED 


Security  Classification 


^CLASSIFIED _ 

Security  Classification 


*CV  WORDS 


Ignition  theory 
Solid  propellant 


INSTRUCTIONS 


ORIGINATING  ACTIVITY:  Enter  thw  name  and  address 
•I  the  contractor,  subcontractor,  gr  imrr,  Department  of  Dr- 
leaae  activity  or  olivet  orfiniiaiua  f  corporate  author)  itaulnf 
the  report. 

2a.  RETORT  SECUOTT  CLASSIFICATION:  Enter  the 
all  security  classification  of  the  report.  Indicate  whether 
’’Restricted  Data**  la  inctodrd  Marking  is  to  be  in  accord* 
•nee  with  appropriate  a  rear  It  y  regulations. 

2b.  CROUP*  A uton at Ic  downgrading  ta  specified  in  DoD  De 
rectlre  5200.10  and  Armed  Forcea  Induatrlal  Manual.  Enter 
the  group  number.  Also,  when  applicable,  ahow  that  optional 
max  king  a  have  been  used  for  Group  J  and  Croup  4  as  author* 
Ued. 

X  REPORT  TITLE  Ertn  the  complete  report  title  in  all 
capital  letters.  Titles  In  all  rasea  should  be  unclassified 
If  a  meaningful  title  cannot  be  selected  without  classifica¬ 
tion,  show  title  classification  In  all  capitals  In  parenthesis 
immediately  following  the  lit.e. 

4.  DESCRIPTIVE  NOTES  If  appropriate,  enter  the  type  of 
report,  e.g.,  interim,  progress,  summary,  annual,  or  final. 

Gtve  the  inclusive  dstes  when  a  specific  reporting  period  ts 
covered 

&.  AUTMORfS):  Enter  the  namefs)  of  autho i<a)  am  shown  on 
or  In  the  report.  Enter  last  name,  first  name,  middle  initial. 

If  nulit ary,  ahow  rank  and  branch  of  service.  The  name  of 
the  principal  author  ia  an  absolute  mi-imum  requirement. 

6.  REPORT  DATEt  Enter  the  dale  of  the  report  as  day. 
month,  year;  or  month,  yean  If  more  than  one  date  appears 
on  the  report,  wse  date  of  publication. 

7a.  TOTAL  NUMBER  OF  PAGES:  The  total  page  count 
should  follow  normal  pagination  procedures,  L e..  enter  the 
number  of  pages  containing  information. 

7k.  NUMBER  OF  REFERENCES  Eater  the  total  number  of 
references  cued  u»  the  rrport. 

•a.  CONTRACT  OR  GRANT  NUMBER:  If  appropriate,  enter 
the  applicable  number  of  the  contract  or  grant  under  which 
the  rrport  was  written 

ftk.  *r.  k  Id  PROJECT  NUMBER:  Enter  the  appropriate 
military  department  Identification,  mach  aa  project  number, 
subprotect  number,  system  numbers,  task  number,  etc. 

9a.  ORIGINATOR’S  REPORT  NUK3ER(SV  Enter  the  offi- 
cial  report  number  by  which  the  document  will  be  identified 
and  controlled  by  the  originating  activity.  This  number  must 
be  unique  to  thiv  report. 

9k.  OTHER  REPORT  NUMBERtS*  If  the  report  has  been 
assigned  any  other  report  numbers  ( either  by  Ike  otigtnmtot 
or  by  ike  sponsor  J.  alao  enter  this  rvuraber(  a„ 

10.  AVAlLAfULfTY/LlMlTATlON  NOTICES:  Enter  any  lim¬ 
itations  on  further  dissemination  of  the  report,  other  than  those 


imposed  by  security  clas'  In  at  ion.  using  standard  statements 
such  as: 

(1)  ”0*alified  requesters  way  obtain  copies  of  this 
report  from  DDC.” 

(2)  ”Fo»  -v  announcement  and  dissemination  of  this 
report  v.  DDC  is  not  authorised.” 

(J)  ”U  S.  Government  agencies  may  obtain  copies  of 
thia  report  directly  from  DDC.  Other  qualified  DDC 
users  shall  request  through 

(4)  ”U.  S  military  agennea  may  obtain  copies  of  this 

report  directly  from  DDC  Other  qualified  users 
shall  request  through 


($)  ”  All  distribution  of  this  report  is  controlled.  ^iaJ- 

ified  DDC  users  shall  request  through 


If  the  repon  has  been  furnished  to  the  Office  of  Technical 
Services.  Depart  me  nl  of  Commerce.  for  sale  in  ihe  public,  indi¬ 
cate  this  fact  and  enter  the  price,  if  known. 

1L  SUPPLEMENTARY  NOTES  Use  for  additional  explana¬ 
tory  notes. 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  ( pay¬ 
ing  tot )  the  research  and  development.  Include  address. 

13  ABSTRACT:  Enter  an  abstract  giving  a  brief  and  factual 
summary  of  the  document  indicative  of  the  rrport.  even  though 
it  r-ay  also  appear  elac where  in  the  body  of  ihe  technical  re- 
port  If  additional  space  is  required,  a  continuation  sheet  shall 
be  attached. 

It  is  highly  desirable  that  the  abstract  of  classified  reports 
be  unclassified.  Each  paragraph  of  the  abstract  shall  end  with 
an  Indication  of  the  military  security  classification  of  the  in¬ 
formation  In  the  paragraph,  represented  as  f  TS).  (Sj.  fCJ.  er  fUj 

There  ts  no  limitation  on  the  length  of  the  abstract.  How¬ 
ever.  the  suggested  length  ia  from  ISO  t?  22Z  words. 

14  KEY  WORD!*.  Key  words  are  technically  meaningful  terms 
or  short  phrases  that  charactenxe  a  report  and  may  be  used  as 
index  entries  for  cataloging  the  report  Key  word*  must  be 
selected  so  that  no  security  classification  is  required.  Identi¬ 
fiers.  such  as  equipment  model  designation,  trade  name,  military 
project  code  name,  geographic  location,  may  be  used  a  a  key 
words  but  anil  be  followed  by  an  indication  of  technical  con¬ 
text.  The  assignment  of  links,  roles,  and  weights  is  optional. 


UNCLASSIFIED 


Security  Classification 


